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Abstract 

In this letter we present a simple relation, that allows one to predict the 
Josephson current between two unconventional superconductors with arbi- 
trary band-structure and pairing symmetry. We illustrate this relation with 
examples of s-wave and d-wave junctions and show a detailed numerical com- 
parison with the phenomenological Sigrist and Rice relation jjj commonly 
used to interpret tunneling experiments between d-wave superconductors. 
Our relation automatically accounts for 'mid-gap states' that occur in d- 
wave superconductors and clearly shows how these states lead to a 

temperature-sensitive contribution to the Josephson current for certain orien- 
tations of the junction. This relation should be useful in exploring the effects 
of disordered junctions, multiple atomic orbitals and complicated order pa- 
rameters 
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I. INTRODUCTION 



Phase-sensitive measurements involving Josephson junctions has recently become an im- 
portant technique for probing the pairing symmetry of high-T c superconductors |5|-0. Sev- 
eral experiments point to the existence of n junctions between high-T c superconductors 
strongly suggesting d-wave symmetry, though other possibilities cannot be ruled out com- 
pletely. In this paper we present a simple expression (valid for junctions with arbitrary 
pairing symmetry), that shows how a particular feature in the density of states contributes 
to the critical current in a Josephson junction. For example, it is now well-known 0] that, for 
certain orientations, d-wave superconductors exhibit a midgap peak in the surface density 
of states. Using our formulation, one can easily see with a 'back of the envelope' calculation 
that this midgap peak will lead to a large contribution to the critical current that will decay 
inversely with temperature. 

We write the current I between the two superconductors 1 and 2 in terms of the phase 
difference <fi as I = I c {4>) sin0, where I c is given by (f(E): Fermi function) : 

Op coo 

I C = ^\M\ 2 dEf(E)J(E) 

(1) 



J{E) =Im[sf<tf e + sftTWi,) 



In this paper we will only discuss the tunneling limit in which I c is independent of and the 
pair correlation functions gf 1 ^ can be approximated by their values at a free surface |§. 
These functions are obtained by evaluating the retarded Green function for the Bogoliubov- 
deGennes equation. Complicated band-structures and pairing symmetries can be accounted 
for by an appropriate choice of the one-particle Hamiltonian H and pair potential A(r, r'). 
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The surface Green function 'g' appearing in eq [I] is obtained from the full Green function 
G(r, r') by setting x = x' = 0, x being the normal to the surface: 

g eh ' h %y,z ]y ',z') = G« h > he {r,r') x=x , =Q 

In eq [I], we have assumed a one-dimensional model for which g eh > he as well as the coupling 
element M are purely numbers. A generalized version of this equation is presented at the 
end of this paper that can be used when g eh > he and M are matrices. For our discussions in 
this paper the simple version in eq[l| is adequate since we will only consider planar junctions 
with s-wave or d-wave order parameters, for which each transverse wave vector (k y , k z ) 
decouples into an independent one-dimensional channel. The generalized version is needed 
to handle disordered junctions, multiple atomic orbitals and complicated multilayer order 
parameters ||||. 

The advantage of this formulation is that it allows us to predict the Josephson current 
from a knowledge of the pair correlation functions for the individual isolated superconduc- 
tors. From this point of view eq [I] is similar to the well-known formula for the conductance 

-^1=4.W^ (3) 

that is routinely used to obtain the density of states (pi^) from experimental data. Indeed, 
the pair correlation functions g eh,he (E) exhibit features similar to those in the density of 
states and can often be guessed intuitively as discussed in section ||. 

We will not go into the derivation of eq [l] which follows from the Green function formal- 
ism ||[10|- However it can be seen easily that for ordinary s-wave superconductors it yields 
the standard results. To establish contact with the literature, we rewrite eq p] as a sum over 
the Matsubara energies uo n = (2n + l)irkT along the imaginary axis: 

h = ^|M| 2 YXatgT + 9?dl\E^ n) (4) 

ft n 

Making use of the property g eh {iuj n ) = g he (iuj n )* = —F(iu n ) (which is true provided 
A(r, r') = A(r', r)), eq § reduces to 



h= ^^\M\ 2 Y,F 2 (iu n )F*{iw n ) (5) 

h n 

This result is often applied to s-wave junctions using the bulk pair correlation function 
for F and can be combined with eq |3| to obtain the Ambegaokar-Baratoff relation 1X1,12 



In d-wave (and some other unconventional) superconductors, the surface pair-correlation 
functions can show features like 'midgap peaks' that are absent in the bulk. The point to 
note is that eq [1| ( or |] ) can be used for arbitrary pairing symmetry, provided we use the 
correct pair- correlation functions at the surface rather than the ones in the bulk. Eqs |1] and 
^ are equivalent and it is largely a matter of taste which one we use. In this paper, we will 
use the former and present our results in terms of the pair correlation functions and the 
current spectrum along the real energy axis. 

II. QUALITATIVE DISCUSSION 

To use eq [l] we need the pair correlation functions g eh . In general it is fairly straight- 
forward to obtain a numerical solution for complicated band-structure and pairing sym- 
metries [pTO | . But even without a detailed quantitative solution, one may be able to guess 



g eh (E) from a knowledge of the density of states and use it to gain insight into the resulting 
effect on the Josephson current. For example, for d-wave superconductors with small values 
of the misorientation angle a and transverse wave- vector (k y ,k z ), g eh has the usual form 
encountered with s-wave superconductors El: 



Jie n eh ^1,2 
01,2 - ^1,2 ~ 



The current spectrum J(E) is s-wave like and the critical current has the usual temperature 
dependence. But when the quantity sin 2a sin k y a is large, g eh is dominated by a singularity 
at E=0 corresponding to the midgap peak in the density of states. 

he eh iAi2 
01,2 = -01,2 \P) 
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This results in a large contribution to the critical current which, however, decays rapidly 
with temperature. To see this, we use eq || in eq |I| to obtain 

J{E) = 2A 1 A 2 -^^ -> 2ixA l A 2 5'{E) as rj -> 

where S'(E) represents the derivative of the delta function. Hence from eq|I| 

T 4 e A A ,^ l2/ df . e A A \M\ 2 



h dE e=o h k B T 

indicating that the contribution of the midgap peak to the critical current decays with 
temperature as l/T. The apparent divergence near T=0 can be eliminated either by 
including a finite rj to represent phase-breaking effects or by going beyond the weak-coupling 
approximation (\Mg eh \ << 1) ||. The recent results of ref |14| for d-wave junctions (which 
includes strong coupling effects) also indicate a similar anomalous temperature dependence 
of the midgap contribution to the critical current. 



III. NUMERICAL RESULTS 

For quantitative calculations, we evaluate the pair-correlation functions from the BdG 
equation using a two-dimensional tight-binding lattice with a dispersion relation of the 
fonn(///Ao = 15, t/A = 10): 

E = —fi + 2i(l — cos k x a) + 2i(l — cosk y a) (7) 

and a pair potential of the form: 

A = Ao cos 2a (cos k x a — cos k y a) 

+A sin(2a)(sin/c a; asinfcj / a) (8) 

We set rj = .025 which corresponds to a phase-relaxation time of = 2 ps, if A is 20 
meV. Fig 2 shows the bulk and surface correlation function for a = and 30 degrees with 
k y a = 0, 1. If either oi or \\iyQj is zero, the midgap term (see eq []) is absent and the surface 
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g looks similar to the bulk g . But with both a and k y a non-zero, the surface and bulk 
functions are completely dissimilar. The bulk pair correlation function exhibits two gaps 
corresponding to +k x and — k x respectively, since they experience different values of A. The 
surface pair correlation function involves a non-trivial mixture of the two and has a midgap 
peak as evident from the second of fig 0b. 

Fig 3a shows the current spectrum for a Josephson spectrum with ai = a 2 = 30 degrees. 
With k y a = 0, the spectrum is s-wave-like with peaks close to E = ±Ai f 2. But with k y a = 1 
we have a large midgap contribution. On integrating over k y ||15| , this midgap contribution 
completely dominates the usual peak away from midgap, (see fig |3]b). The midgap peak 
has the form of a doublet function (5'(E)) and its contribution to I c will disappear with 
increasing temperature as discussed earlier. 

Fig [|a and b show the critical current as a function of the misalignment angles «i and 
«2 with kT = and with kT = .5 A respectively. For comparison we have also shown the 
results from the phenomenological relation given by Sigrist and Rice jl] (Fig [|c): 

I c — Ico cos 2«i cos 2a 2 (9) 

We note that in fig f|a, the current has peaks at a>\ = ±7r/4,a;2 = ±vr/4 due to the large 
midgap peaks in g eh , in sharp disagreement with eq |9] |16|]. At higher temperatures (fig 
|]b) the contribution of these peaks is reduced as explained earlier, and the variation in the 
critical current looks closer to what we expect from eq|9|. A full calculation of I C {T) requires 
a self-consistent evaluation of the pair potential ||, which we do not address in this paper. 
Using the self-consistent pair potential from ref ||, we find that the current contribution 
due to the midgap peak is suppressed somewhat, but its effect should still be observable. 

IV. GENERALIZATION 

Finally we note that eq [1| can be generalized to include cases where the pair correlation 
functions as well as the matrix element M are matrices instead of pure numbers. The current 
spectrum is then given by 
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J{E) = lm[Tr[g e 2 h M + g^M + g^M + g\ h M}} {E+lv) (10) 

instead of eq|l|. We believe eq[l(] should be useful in calculating the Josephson current taking 
into account disordered junctions, multiple atomic orbitals, or complicated pair potentials. 
The required surface Green functions can be calculated numerically using an iterative scheme 
once the appropriate tight-binding Hamiltonian has been identified [|l(J . 
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FIGURES 
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FIG. 1. Structure considered in this paper. M is the matrix element linking superconductors 
1 and 2. In case of d-wave, the pair potential in each superconductor varies with the direction 
9 according to A(6>) = Aocos(2# — 2a), where a is the angle that the main lobe makes with the 
x-axis. 
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(b) a=30 
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FIG. 2. Magnitude and phase of the pair correlation function g eh at a surface for a = degrees 
and 30 degrees with (i) k y a = and (ii) k y a = 1. These correspond to (i) 6* = and (ii) = 50 
degrees respectively, since sin# = k y /kf and fe/a = 1.3 . Also shown is the magnitude of the bulk 
pair correlation function. ( — ^surface(magnitude), : surface (phase), xxx: bulk(magnitude) ) 
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(b) integrated over all k y a 



A J(E) 
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FIG. 3. Josephson current spectrum for a d-wave junction having ol\ = c*2 = 30 degrees (a) 
with k y a = 0, 1 (xxx : k y a = 0, solid : k y a = 1) , (b) integrated over all values of k y a. 
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FIG. 4. Critical current (normalized) of Josephson junctions as a function of the misalignment 



angles ot\ and 0:2 (in units of it): (a) kT = 0, (b) kT = .5Aq and (c) Phenomenological result based 



on I c = I C QCOs(2ai)cos(2a2) 
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